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Abstract. - We study the probability, Ps{t), of a cluster to remain intact in one-dimensional 
cluster-cluster aggregation when the cluster diftusion coefficient scales with size as D{s) ~ s^. 
Ps(t) exhibits a stretched exponential decay for 7 < and the power-laws t~^^^ for 7 = 0, and 
^-2/(2 I-) £qj. < 7 < 2. A random walk picture explains the discontinuous and non-monotonic 
behavior of the exponent. The decay of Ps{t) determines the polydispersity exponent, r, 
which describes the size distribution for small clusters. Surprisingly, t{"/) is a constant r = 
for < 7 < 2. 



Many models of aggregation phenomena lead to scale-invariance: the average cluster size 
increases as a power-law, S{t) ~ t^, which defines a dynamical exponent z. This kind of 
behavior is met in various contexts ranging from chemical engineering to materials science 
to atmosphere research to, ultimately, even astrophysics It is of interest to explore the 
statistics of aggregation as a dynamical process, beyond the length- and timescales defined 
through z. 

In this Letter we introduce a new quantity in aggregation systems, the cluster survival, 
defined as the probability Ps{t), that a cluster present at t = remains unaggregated until 
time t. This is a first passage problem j2| in a many body system and analogous to persistence 
which is often studied by measuring the fraction of a system that preserves its initial condition 
for all times [0,t] The cluster survival turns out to decay in a nontrivial and counterin- 
tuitive manner. The behavior can be understood by a mean-field like random walk analysis. 
However, even on the mean-field level the question reduces to a novel, unsolved random walk 
problem, which we analyse in the long time limit. More importantly, by solving the decay of 
the cluster survival we are able to determine the polydispersity exponent characterising the 
cluster size distribution. 

We concentrate on a common and important example: diffusion-limited cluster-cluster 
aggregation (DLCA) Q . In the lattice version of DLC A any set of nearest neighbor occupied 
lattice sites is identified as a cluster. Each of these performs a random walk with a size 
dependent diffusion constant, D{s) ~ s'^, where 7 is the diffusion exponent. Colliding clusters 
are merged together and the aggregate diffuses either faster (7 > 0) or slower than before 
(7 < 0). In the following cluster survival is investigated in the one-dimensional case for 
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numerical and analytical simplicity. We employ numerical simulations and random walk (RW) 
arguments. 

The numerics is made transparent by mapping the behavior of surviving clusters to a 
three-particle RW picture: two particles with a time-dependent diffusion coefficient confine a 
surviving one which diffuses at a constant rate. This is an analogy of the famous independent 
interval approximation often used in persistence studies 1^. One can discern three separate 
cases: first, < 7 < 2, which results in a power-law decay for the survival: Psit) ~ t~^s(7)^ 
second, 7 = which is exactly solvable both for the RW and the DLCA problems, and 7 < 
when Ps decays as a stretched exponential. For 7 > 2 the system has a gelation transition 
and is not of interest here. 

The RW problem has, to our knowledge, not been discussed in the literature. We consider 
the asymptotic behavior of the associated Fokker-Planck equation and perform numerical 
simulations. The first main result is that 6*3(7) — 2/(2 — 7) when < 7 < 2. Thus the 
survival is discontinuous and non-monotonic since 9s{0) =3/2 but ^5(0+) = f . This non- 
intuitive result follows since for 7 > the RW problem becomes asymptotically separable 
as the ratio of the diffusion coefficients diverges: the surviving clusters are immobile. For 
7 < the same is not true and the fluctuations of the confining clusters remain relevant in 
determining the stretching exponent. 

One of the main interests in aggregation is the behavior of the cluster size distribution, 
ns{t) (the number of cluster of size s per lattice site at time t). For DLCA simulations and 
experiments have validated the scaling Q ns{t) = S{t)~^f {s/S{t)) , where the scaling limit, 
s — > 00, S{t) — > 00 with X = s/S{t) fixed, is taken. In one dimension z = 1/(2 — 7) [^. 

There is a fundamental difference between 7 < and 7 > 0. For 7 < the cluster 
size distribution is bell-shaped and f{x) decays faster than any power at both tails. For 
7 > it is broad so that f{x) ~ x^'^ {x 0) defines the polydispersity exponent, t, which 
characterizes the density of small clusters. In this region we show evidence for the scaling 
relation 0s = {2 — t)z. Therefore r is determined by the cluster survival strategy, shedding 
light on the non-trivial problem how to compute it The final main result follows thus: 

t(7) — for < 7 < 2 indicating a flat cluster size distribution. This is conflrmed by 
simulations. 

Next we present a mean-field random walk analysis to calculate the survival exponent 
9s{j)- Although the kinetics in one-dimension is fluctuation-dominated the mean-field 
approach turns out to capture the essential ingredients for 7 > 0. This is demonstrated by 
comparing the random walk survival to that of the full DLCA one. In the opposite case, when 
7 < 0, the stretching exponent obtained from RW simulations differs from the DLCA case. In 
the limit t 3> the average distance between clusters grows as and between the surviving 
ones as t^^ . The latter become separated by aggregated clusters at late times, since evidently 
6s > z. Thus it is sufficient to consider only one trial cluster and its two neighbors. These 
grow still by collisions with neighbors on the opposite side. In the mean-field approximation 
the discrete growth events can be substituted by a continuous process and the neighboring 
clusters grow as the average cluster does. The finite extent of clusters is irrelevant and they 
can be considered as point particles. Let Xi{t) [i — 1,2,3) denote their positions at time t 
with a;i(0) < 2:2(0) < 2:3(0). 

The motion of these particles is described by the Langevin equations 

i^{t)^m (1) 

with Gaussian white noises {£,i{t)) = and (Ci(^)Ci(^')) = 2'Di{t)SijS{t — t'), in the standard 
notation. The diffusion coefficients read as Vi{t) = 2?3(t) = Dif^^ and T>2{t) = 1)2 ■ The time 
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dependent diffusion constant, say I'i(t), implies that the particle f will follow a simple diffusive 
motion with a constant diffusion coefficient Di in a time scale T{t) — dt''Di{t')/ Di = 

The survival of the center particle is determined by the termination of the process, given 
by either xi{t) = X2{t) or 2:2 (i) = 2^3 (i)- It is natural to consider the distances between the 
particles: Xi2{t) — X2{t) — xi{t) > and X23{t) — x^it) — X2{t) > 0. Starting from Eq. |^ the 
following Fokker-Planck equation is reached: 



dt \9xl2 dx^^J 8x128x23' 

where p{xi2, X23'-,t) is the probability density for the two distances at time t. The initial 
condition is now p{xi2, X23'-,0) = 5{xi2 — a;i2(0))(S(a;23 — 2^23(0))- The termination of the 
process when two particles collide gives absorbing boundary conditions along the axis, i.e., 
p{xi2, 0; = and p(0, X23] t) = Q for all times t. The survival probability 

/•oo /"OC 

Pnwit)^ dxi2 dx23 p{xi2,X23;t)^t-''^^'^'^\ (3) 







where the exponent 9rw is the survival exponent for the RW problem. 

Equation |^ can not be solved for t arbitrary since the absorbing boundary conditions to- 
gether with the two different time scales make the standard methods inappropriate. However, 
the survival exponent is given by the leading-order asymptotic behavior when i — > oo. We 
consider the large time limit and the three different cases separately: 7 < 0, 7 = 0, and 7 > 0. 
In the size independent case, 7 = 0, the collisions of the clusters surrounding a surviving 
cluster with other clusters do not matter. This is an old problem of three similar annihilating 
random walkers, for which the survival exponent is known to be ORwi^) = 3/2 |lOj. Also the 
full DLCA problem can be solved exactly with the result 0s{O) — 3/2 pl] |. 

Diagonalizing Eq. ^ in the time scale T gives a simple diffusion equation in a wedge with 
two absorbing boundaries. However, for 7^0 the wedge angle is a function of time making 
the exact solution hard. For 7 > the leading term of the survival probability at late times 
can be obtained by considering the survival in the final wedge angle Ooo = 7^/2. Thus Q 
-Prw(0 ^ T-''/26'°o ^ T-^ - Since z = 1/(2 - 7) the survival exponent reads 

6rw{i) = 2z = 2/(2 — 7). The approximation obtained by replacing the time-dependent 
angle by the final opening angle corresponds to putting D2 = va Eq. H, i.e., to taking the 
center particle to be at rest. This guess can be validated by directly solving equation g and 
analyzing the limiting behavior of the solution Note, that in this limit 9Ry/{l) can be 
simply determined from two independent random walkers with a fixed absorbing boundary in 
between. 

For 7 < the situation is more tricky and will be considered in more detail elsewhere ||l^ . 
Briefly, proceeding similarly as above leads to a closing wedge with final angle 800 = 0. This 
would correspond to a situation where now the particles 1 and 3 are fixed and therefore to a 
simple exponential decay for the survival ||^. However, simulations show that the expectation 
value of their positions grows like t"' with non-trivial a. The survival becomes a stretched 
exponential PRw(i) ~ exp(— Ci''^™) with C > a constant. If the distance between the out- 
ermost particles would grow deterministically as t" this would lead to a stretched exponential 
decay with /3det = 1 — 2a Here the fluctuations in the particle position violate this relation. 

The RW picture is tested by comparing it to DLCA by simulations. These are done in the 
usual fashion, with a monodisperse initial condition and equal distances between neighboring 
clusters. Pg is independent of the initial distribution except for transient effects [Q. 
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In the case 7 = with the initial distances being xi2(0) = 2:23(0) = the first correction 

3/2 



O 



'0 

Dt 



SO that the cor- 



to scahng reads Ps{t) = ^ (4) x |l - ^ 

rection becomes neghgible for times much larger than the cross-over time ta — 31q/ {16Dt). 
For 7 > the ratio of the diffusion coefhcients, r = D2/Dit'^^, controls the validity of the ap- 
proximation of neglecting the constant terms in Eq. |^. Therefore the cross-over time depends 
on 7 as tcr ^ r'^'^-'i)h ^ which diverges for 7^0. We can thus expect that the asymptotic 
scaling regime can be reached in simulations only for relatively large values of 7. 

Figure |l| shows the survival probabilities obtained from simulations. Prw(^) clearly decays 
as a power-law for large times. For 7 = the survival exponent saturates to the asymptotic 
value Ojiw — 3/2 around t k, 10^ as shown in the inset, where the local survival exponent is 
presented. The exponent saturates also for 7 — 0.75 and 7 = 1. For < 7 < 0.75 the survival 
exponent is slowly approaching the asymptotic value given by Orw — 2/(2 — 7). For 7 = 0.75 
r « 60 at the time when the exponent saturates. This would correspond to icr ~ 2 x 10^ and 
3 X 10^^ for 7 = 0.50 and 0.25, respectively. Thus we can not reach the asymptotic regime for 
7 ^ 0.5. 

In Figure |^ the DLCA and RW survival probabilities are compared. The two behaviors 
agree except for a small difference between the amplitudes. The initial inter-particle distances 
are taken to be the same, in order the RW-picture to be as close to DLCA as possible. Notice 
the figure contains a case (7 = 0.33) in which the asymptotic regime is not reached. 

For 7 < the survival probability decays stretched exponentially for both RW and 
DLCA systems (not shown). However, the stretching exponents differ from each other. 
This further supports the importance of fluctuations for 7 < 0. For RW survival we ob- 
tain /3rw ~ 0.33,0.31,0.25, and 0.19 for 7 = -6,-4,-2, and -1, respectively. For DLCA 
the numerics suggest an expression Ps = 2/3(1 — 2z). We also separately checked that the 
average size of the neighbors of unaggregated clusters in DLCA grows as . Note, that 
lim-),_>_oo Ps = 2/3 and not a simple exponential {j3s = 1) as one might expect based on the 
case of two immobile neighbors. It is also quite surprising, that the mean field approximation 
works better for a broad cluster size distribution (7 > : f{x) ^ x^'^{x 0)) than in the 
case, where the distribution is narrow around the mean (7 < ; f{x) decays faster than any 
power law at both ends). 

To summarize the results the survival probability decays as 



Ps{t)-{t-^/' ,7-0 (4) 




The survival exponent is discontinuous at 7 = 0, i.e., 3/2 — 9s{0) > Os{0^) — 1. This seems 
first counterintuitive since making some of the particles to diffuse faster helps the other to 
survive longer! This is the best since the surviving particle eventually "discovers" the optimal 
strategy of remaining stationary. This is further confirmed by the fact that the probability 
of finding a site that has never been under any of the clusters decays as power-law with the 
exponent 2/(2 — 7) for all 7 < 2 [|l3|. For 7 < the surviving clusters no longer remain 
stationary. 

The third exponent of interest in DLCA scaling is the decay exponent, w, which describes 
the decrease of the number of clusters of a fixed size s as a function of time Us (t) ~ . As 
the other exponents z and r, it is expected to depend on 7. The three exponents are related 
by the scaling relation |l6[ w = (2 — t)z. Therefore dynamic scaling is fully characterized by 
any two of the exponents. Even on mean-field level (Smoluchowski's equation) the only easy 



E. Hellen: Particle survival 



5 



exponent is z since it does not involve the full scaling function as in the case of and t and 
w [|. 

For 7 = an exact solution of the cluster size distribution na(t) is possible, with w = 3/2 for 
any short-range correlated initial distribution (0) JlTf . For a monodisperse initial condition, 
ns(0) = (5i,s, the survival probability is simply ni(t), yielding the exponent 0s(O) = 3/2. For 
7 > the dynamics of the clusters in the s <^ S{t) part of the size distribution is dictated 
by collisions with larger, faster ones, that remove such small clusters from the tail. The 
mechanism by which clusters stay in the tail should be the same as for the survival problem. 
Thus for < 7 < 2 one should have for the decay exponent w — 9 s = 2 / {2 — — 2z together 
with the scaling relation 

es = {2- t)z. (5) 

The numerically estimated values for the exponents fulfill Eq. ^ within the error bars for all 
values of the diffusion exponent 7 . Hence, the polydispersity exponent is discontinuous at 
7 = since t(0) — —1 ^ Q = t(0"'"), similarly to some examples on the mean- field level p7| . 
It is also surprisingly enough independent of the value of 7. 

Simulations confirm this although crossover effects make the analysis intractable near 
7 = 0. Figure ^ shows the scaling plots for the cluster size distribution for various values 
of the diffusion exponent. The bigger the 7 is the faster the scaling function approaches a 
constant near x — Q. For 7 < 0.5 the times reached in simulations are too small to reveal the 
asymptotic scaling behavior. In this region also the measurements of w from 7is(t) ^ t^^ give 
information only on the crossover effects . 

In summary, we have studied the survival of clusters in DLCA in one dimension. The 
decay of the initial state equals the density of clusters that stay intact by not aggregating 
with others. This can be analyzed on the mean-field level as the survival of a random walker 
bounded by two others with time-dependent diffusion coefficients. This maps to diffusion in a 
wedge with absorbing boundaries and time-dependent wedge angle. For 7 > the surviving 
particles are such that they pick the strategy of staying immobile. The resulting survival 
exponent is non-monotonic and discontinuous at 7 = 0. Cluster survival determines the 
polydispersity exponent r that characterizes the small cluster size tail. It is discontinuous 
and constant, r = 0, for < 7 < 2. For 7 < the survival probability decays stretched 
exponentially, the fluctuations of the neighboring clusters determine the stretching exponent, 
and the mean-field RW-picture gives only a qualitative understanding of the survival. 

Above one dimension the study of cluster survival will be both interesting and much less 
straightforward. A similar mean-field picture in terms of first-passage times of random walks 
is not directly applicable. For example, one lacks much of the theory needed to analyze 
the survival behavior of many interacting particles. We conclude with the conjecture that 
the solution of the survival problem is related to the cluster size distribution also in higher 
dimensions. There it should be possible to find experimental realizations, to study the survival 
phenomenon p^ . 

Acknowledgements - We thank Erik Aurell and Amit K. Chattopadhyay for discussions. 
This research has been supported by the Academy of Finland's Center of Excellence program. 

REFERENCES 

[1] See, for example. Kinetics of Aggregation and Gelation, edited by F. Family and D. P. Landau 

(North-Holland, Amsterdam, 1984). 
[2] S. Redner, A Guide to First-Passage Processes (Cambridge University Press, New York, 2001). 



6 



EUROPHYSICS LETTERS 



10 k/l/lK' 



10" 



tr 



10" 



10" 





0.00 




0.25 




0.50 




0.75 


- Y = 


1.00 



10" 




Fig. 1 - The survival probability of three random walkers. The inset shows the corresponding local 
exponents. The solid lines correspond to the analytic values given by 6rw = 2/(2 — 7). The data are 
averaged over variable number of realizations ranging from 10^ for 7 = to 2 X 10^ for 7 = 0.5. The 
initial distance between particles is 10. 



[3] S. N. Majumdar, Curr. Sci. (India) 77, 370 (1999). 

[4] P. Meakin, Phys. Scripta 46, 295 (1992). 

[5] S. N. Majumdar, C. Sire, A. J. Bray, and S. J. Cornell, Phys. Rev. Lett. 77, 2867 (1996); B. 
Derrida, V. Hakim, and R. Zeitak, Phys. Rev. Lett. 77, 2871 (1996). 

[6] K. Kang, S. Redner, P. Meakin, and F. Leyvraz, Phys. Rev. A 33, 1171 (1986); S. Miyazima, 
P. Meakin, and F. Family, Phys. Rev. A 36, 1421 (1987). 

[7] P. G. J. van Dongen and M. H. Ernst, Phys. Rev. Lett. 54, 1396 (1985). 

[8] S. Cueille and C. Sire, Phys. Rev. E 55, 5465 (1997). 

[9] K. Kang and S. Redner, Phys. Rev. A 30, 2833 (1984). 

[10] M. E. Fisher, J. Stat. Phys. 34, 667 (1984). 

[11] J. L. Sponge, Phys. Rev. Lett. 60, 871 (1988). 

[12] E. K. O. Hellen, P. E. Salmi, and M. J. Alava, in preparation. 

[13] E. K. O. Hellen and M. J. Alava, in preparation. 

[14] B. Derrida and R. Zeitak, Phys. Rev. E 54, 2513 (1996). 

[15] S. Redner and P. L. Krapivsky, Am. J. Phys 67, 1277 (1999). 

[16] T. Vicsek and F. Family, Phys. Rev. Lett. 52, 1669 (1984). 

[17] P. G. J. van Dongen and M. H. Ernst, Phys. Rev. A 32, 670 (1985). 

[18] E. K. O. Hellen, T. P. Simula, and M. J. Alava, Phys. Rev. E 62, 4752 (2000). 



E. Hellen: Particle survival 



7 




7=0.33; RW 
□O 7=0.33; DLCA 
RW 
DLCA 



10^ 



10^ 



10^ 



10^ 



Fig. 2 - a) Comparison between the random walk and DLCA survival probabilities, b) The corre- 
sponding local exponents. The solid lines correspond to the analytic values given by 2/(2 — 7). The 
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The DLCA simulations are averaged over 50000 realizations 
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Fig. 3 - The scaling of the cluster size distributions for 7 = 0.00, 0.40, 0,57, and 1,00 (from bottom 
to top). Distributions are shown at times from f = 2^ to 2'°'"'"', where fcmax = 17, 15, 14, and 11 for 
7 = 0.00, 0.40, 0,57, and 1,00, respectively. The dashed line has a slope 1. The data are averaged 
over 50000 realizations on a system of size 55555. The distributions (except the 7 = distribution) 
have been shifted in the vertical direction for clarity. 



